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q ■ Let P be a set of n points in M 3 . The 2-center problem for P is to find two congruent 

balls of minimum radius whose union covers P. We present two randomized algorithms 
for computing a 2-center of P. The first algorithm runs in 0(n 3 log 5 n) expected time, 
and the second algorithm runs in 0((n 2 log 5 n)/(l — r*/ro) 3 ) expected time, where r* 
is the radius of the 2-center balls of P and ro is the radius of the smallest enclosing 
ball of P. The second algorithm is faster than the first one as long as r* is not too 
close to ro, which is equivalent to the condition that the centers of the two covering 
balls be not too close to each other. 
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1 Introduction 



1.1 Background 

Let P = {pi, . . . ,p n } be a set of n points in M 3 . The 2-center problem for P is to find 
two congruent balls of minimum radius whose union covers P. This is a special case of 
the general p-center problem in M. d , which calls for covering a set P of n points in M. d by p 
congruent balls of minimum radius. If p is part of the input, the problem is known to be 
NP-complete [29] even for d = 2, so the complexity of algorithms for solving the p-center 
problem, for any fixed p, is expected to increase more than polynomially in p. Agarwal 
and Procopiuc showed that the p-center problem in M. d can be solved in n°( pl time [2], 
improving upon a naive solution. At the other extreme end, the 1-center problem (also 
known as the smallest enclosing ball problem) is known to be an LP-Type problem, and 
can thus be solved in 0(n) randomized expected time in any fixed dimension, and also in 
deterministic linear time [T^J [27J I2B]- Faster approximate solutions to the general p-center 
problem have also been proposed [21 [HE]. 

If d is not fixed, the 2-center problem in R d is NP- Complete [30]. The 2-center problem 
in IR 2 has a relatively rich history, mostly in the past two decades. Hershberger and Suri [23] 
showed that the decision problem of determining whether P can be covered by two disks of a 
given radius r can be solved in 0(n 2 logn) time. This has led to several nearly-quadratic al- 
gorithms [31 [201 123] that solve the optimization problem, the best of which, due to Jaromczyk 
and Kowaluk [23], runs in 0(n 2 log n) deterministic time. Sharir [34] considerably improved 
these bounds and obtained a deterministic algorithm with 0(nlog 9 n) running time. His al- 
gorithm combines several geometric techniques, including parametric searching, searching in 
monotone matrices, and dynamic maintenance of planar configurations. Chan [T2] (following 
an improvement by Eppstein [21]) improved the running time to O (n log 2 n log 2 logn). 

The only earlier work on the 2-center problem in IR 3 we are aware of is by Agarwal et al. [1] , 
which presents an algorithm with 0(n 3+£ ) running time, for any e > 0. It uses a rather com- 
plicated data structure for dynamically maintaining upper and lower envelopes of bivariate 
functions. 

1.2 Our results 

We present two randomized algorithms for the 2-center problem in R 3 . We first present 
an algorithm whose expected running time is 0(n 3 log 5 n). It is conceptually a natu- 
ral generalization of the earlier algorithms for the planar 2-center problem 0, [201 123]; 
its implementation however is considerably more involved. The second algorithm runs in 
0((n 2 log 5 n)/(l — r*/Vo) 3 ) expected time, where r* is the common radius of the 2-center 
balls and r$ is the radius of the smallest enclosing ball of P. This is based on some of the 
ideas in Sharir's planar algorithm [33], but requires several new techniques. As in the previ- 
ous algorithms, we first present algorithms for the decision problem: given r > 0, determine 
whether P can be covered by two balls of radius r. We then combine it with an adaptation 
of Chan's randomized optimization technique [H] to obtain a solution for the optimization 
problem. In both cases, the asymptotic expected running time of the optimization algorithm 
is the same as that of the decision procedure (which itself is deterministic). 
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The paper is organized as follows. Section [2] briefly sketches our two solutions. Section [3] 
presents the near-cubic algorithm, and Section H] presents the improved algorithm. A key 
ingredient of both algorithms is a dynamic procedure for testing whether the intersection of 
a collection of balls in M 3 is nonempty. We present the somewhat technical details of this 
procedure in Section |5l and conclude in Section [6] with a few open problems. 

2 Sketches of the Solutions 
2.1 The near-cubic algorithm 

To solve the decision problem, in the less efficient but conceptually simpler manner, we 
use a standard point-plane duality, and replace each point p G P by a dual plane p*, and 
each plane h by a dual point h*, such that the above-below relations between points and 
planes are preserved. We note that if P can be covered by two balls Bi, B 2 (not necessarily 
congruent), then there exists a plane h (containing the circle dBi (1 dB 2 , if they intersect 
at all, or separating Bi and B 2 otherwise) separating P into two subsets Pi,P 2 , such that 
P\ C B\ and P 2 C B 2 . We therefore construct the arrangement A of the set {p* \ p G P} of 
dual planes. It has 0(n 3 ) cells, and each cell r has the property that, for any point w G r, its 
primal plane w* separates P into two subsets of points, P+ and P~, which are the same for 
every w G r, and depend only on r. We thus perform a traversal of A, which proceeds from 
each visited cell to a neighbor cell. When we visit a cell r, we check whether the subsets 
P+ and P~ can be covered by two balls of radius r, respectively. To do so, we maintain 
dynamically the intersection of the sets {B r (p) \ p G P^ }, {B r (p) \ p G P~}, where B r (p) is 
the ball of radius r centered at p, and observe that (a) any point in the first (resp., second) 
intersection can serve as the center of a ball of radius r which contains P+ (resp., P~), and 
(b) no ball of radius r can cover P+ (resp., P~) if the corresponding intersection is empty. 
Moreover, when we cross from a cell r to a neighbor cell r', P+ changes by the insertion or 
deletion of a single point, and P~ undergoes the opposite change, so each of the sets of balls 
{B r {p) | p G P+ }, {B r {p) | p G P~} changes by the deletion or insertion of a single ball. As 
we know the sequence of updates in advance, maintaining dynamically the intersection of 
either of these sets of balls can be done in an offline manner. Still, the actual implementation 
is fairly complicated. It is performed using a variant of the mult i- dimensional parametric 
searching technique of Matousek [26] (see also [TQl HZl El] ) . The same procedure is also used 
by the second improved algorithm. For the sake of readability, we describe this procedure 
towards the end of the paper, in Section 

The main algorithm uses a segment tree to represent the sets P+ (and another segment 
tree for the sets P~). Roughly, viewing the traversal of A as a sequence £ of cells, each ball 
B r (p) has a life-span (in P+), which is a union of contiguous maximal subsequences of cells r, 
in which p G P+, and a complementary life-span in P~ . We store these (connected portions 
of the) life-spans as segments in the segment tree. Each leaf of the tree represents a cell r of 
A, and the balls stored at the nodes on the path to the leaf from the root are exactly those 
whose centers belong to the set P+ (or P~). By precomputing the intersection of the balls 
stored at each node of the tree, we can express each of the intersections f]{B r (p) | p G P+} 
and f]{B r (p) | p G P~}, for each cell r, as the intersection of a logarithmic number of 



2 



precomputed intersections (see also [20]). We show that such an intersection can be tested 
for emptiness in 0(log 5 n) time. This in turn allows us to execute the decision procedure 
with a total cost of 0(n 3 log 5 n). We then return to the original optimization problem and 
apply a variant of Chan's randomization technique [11] to solve the optimization problem 
by a small number of calls to the decision problem, obtaining an overall algorithm with 
0(n 3 log 5 n) expected running timeQ 

2.2 The improved solution 

The above algorithm runs in nearly cubic time because it has to traverse the entire arrange- 
ment A, whose complexity is 0(n 3 ). In Section H] we improve this bound by traversing only 
portions of A, adapting some of the ideas in Sharir's improved solution for the planar prob- 
lem [33]. Specifically, Sharir's algorithm solves the decision problem (for a given radius r) in 
three steps, treating separately three subcases, in which the centers ci, C2 of the two covering 
balls are, respectively, far apart ( | ci c 2 1 > 3r), at medium distance apart (r < \c\C2\ < 3r) 
and near each other (|ciC2| < r). We base our solution on the techniques used in the first 
two cases, which, for simplicity, we merge into a single case (as done in [2T] for the planar 
case), and extend it so that we only need to assume that \c\C2\ > /3r, for any fixed (3 > 0. 
In more detail, letting B r (p) denote the disk of radius r centered at a point p, Sharir's algo- 
rithm guesses a constant number of lines I, one of which separates the centers c±,C2 of the 
respective solution disks Di,D2, so that the set Pl of the points to the left of I is contained 
in D\. We then compute the intersection K(P L ) = f] p ep L B r (p), and intersect each dB r (p), 
for p G Pr = P \ Pl (the subset of points to the right of /), with dK(P L ). It is easily seen 
that dK(Pi) has linear complexity and that each circle dB r (p), for p e Pr, intersects it at 
two points (at most). This produces 0(n) critical points (vertices and intersection points) 
on 8K(Pl) and 0(n) arcs in between. As argued in [33], it suffices to search these points 
and arcs for possible locations of the center of Di (and dynamically test whether the balls 
centered at the uncovered points have nonempty intersection). 

Generalizing this approach to 1R 3 , we need to guess a separating plane A, to retrieve the 
subset Pl C P of points to the left of A, to compute OK^Pl) (which, fortunately, still has 
only linear complexity), to intersect dB r (p), for each p e Pr, with OK(Pl), and to form the 
arrangement of the resulting intersection curves. Each cell of this arrangement is a candidate 
for the location of the center of the left covering ball B 1 , and for each placement in r, Bi 
contains the same fixed subset of P (which depends only on r). 

However, the complexity of the resulting arrangement Mk on 8K(Pl) might potentially 
be cubic. We therefore compute only a portion M of Mk, which suffices for our purposes, 
and prove that its complexity is only 0(n 2 ). This is the main geometric insight in the 
improved algorithm, and is highlighted in Lemma 14.11 We show that if there is a solution 
then 0(1/ f3 3 ) guesses suffice to find a separating plane. This implies that the running time 
of the improved decision procedure is 0((l//3 3 )n 2 log 5 n). Thus, it is nearly quadratic for 
any fixed value of j3. We show that one can take (3 = 2(ro/r — 1), where tq is the radius of 

1 The earlier algorithm in [1] follows the same general approach, but uses an even more com- 
plicated, and slightly less efficient machinery for dynamic emptiness testing of the intersection of 
congruent balls. 
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the smallest enclosing ball of P. 

To solve the optimization problem, we conduct a search on the optimal radius r*, using 
our decision procedure, starting from small values of r and going up, halving the gap between 
r and tq at each step[§, until the first time we reach a value r > r*. Then we use a variant 
of Chan's technique [11], combined with our decision procedure, to find the exact value of 
r* . The way the search is conducted guarantees that its cost does not exceed the bound 
0((1/ (3 3 )n 2 log 5 n), for the separation parameter (3 = 2(r /r* — 1) for r*. Hence, we obtain a 
randomized algorithm that solves the 2-center problem for any positive separation of C\ and 
C2, and runs in 0((n 2 log 5 n)/(l — r*/r ) 3 ) expected time. 

3 A Nearly Cubic Algorithm 

3.1 The decision procedure 

In this section we give details of the implementation of our less efficient solution, some of 
which are also applicable for the improved solution. Recall from the description in Section [2] 
that the decision procedure, on a given radius r, constructs two segment trees T + ,T~, on 
the life-spans of the balls B r (p), for p £ P (with respect to the tour of the dual plane 
arrangement A). Each leaf is a cell r of A, and the balls, whose centers belong to (resp., 
P~), are those stored at nodes on the path from the root to r in T + (resp., T~). 

For each node u of T + , let S u denote the intersection of all the balls (of radius r) stored 
at u. We refer to each S u as a spherical polytope; see [6j [7J [8] for (unrelated) studies of 
spherical polytopes. We compute each S u in 0(|*S' lt | log \S U \) deterministic time, using the 
algorithm by Bronnimann et al. [5] (see also [HI E2] for alternative algorithms). Since the 
arrangement A consists of 0(n 3 ) cells, standard properties of segment trees imply that the 
two trees require 0(n 3 logn) storage and 0(n 3 log 2 n) preproccessing time. 

Clearly, the intersection K(P+) (resp., K(P~)) of the balls whose centers belong to 
(resp., P~) is the intersection of all the spherical polytopes S u , over the nodes u on the path 
from the root to r in T + (resp., T~). 

Intersection of spherical polytopes. Let S = {Si, . . . , St} be the set of t = O(logn) 
spherical polytopes stored at the nodes of a path from the root to a leaf of T + or of T~, 
where, as above, a spherical polytope is the intersection of a finite set of balls, all having 
the common radius r. Each Si is the intersection of some rij balls, and ^21=1 n i — n - Our 
current goal is to determine, in polylogarithmic time, whether the intersection K of the 
spherical polytopes in S is nonempty. If this is the case for at least one path of T + and for 
the same path in T~ then r* < r, and otherwise r* > r. Moreover, if there exist a pair of 
such paths for which both intersections have nonempty interior, then r* < r (because we can 
then slightly shrink the balls and still get a nonempty intersection). If no such pair of paths 
have this property, but there exist pairs with nonempty intersections (with at least one of 
them being degenerate) then r* = r. 

2 We have to act in this manner to make sure that we do not call the decision procedure with 
values of r which are too close to ro, thereby losing control over the running time. 
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The algorithm for testing emptiness of K is technical and fairly involved. For the sake 
of readability, we delegate its description to Section [5J It uses a variant of multidimensional 
parametric searching which somewhat resembles similar techniques used in PJJJ [HJ HHJ ETJ • 
It is essentially independent of the rest of the algorithm (with some exceptions, noted later). 
We summarize it in the following proposition. 

Proposition 3.1. Let S be a collection of spherical polytopes, each defined as the intersection 
of at most n balls of a fixed radius r. Let N denote the sum, over the polytopes of S, of the 
number of balls defining each polytope. After a preprocessing stage, which takes O(Nlogn) 
time and uses O(N) storage, we can test whether any t < logn polytopes of S have a 
nonempty intersection in 0(log 5 n) time, and also determine whether the intersection has 
nonempty interior. 

Hence, we check, for each cell r, whether each of K(P+) and K(P~) are nonempty and 
non-degenerate. To this end, we go over each path of T + , and over the same path of T~, 
and check, using the procedure described in Proposition [3J], whether the spherical polytopes 
along the tested paths (of T + and of T~) have a nonempty intersection (and whether these 
intersections have nonempty interiors). We stop when a solution for which both K(P+) and 
K(P~) are nonempty and non-degenerate is obtained, and report that r* < r. Otherwise, 
we continue to test all cells r. If at least one degenerate solution is found (i.e., a solution 
where both K(P+), K(P~) are nonempty, and at least one of them has nonempty interior), 
we report that r* = r, and otherwise r* > r. 

By proposition 13. 1[ the cost of this procedure is 0(n 3 log 5 n). This subsumes the cost of 
all the other steps, such as constructing the arrangement A and the segment trees T + ,T~. 
We therefore get a decision procedure which runs in 0(n 3 log 5 n) (deterministic) time. 

3.2 Solving the optimization problem 

We now combine our decision procedure with the randomized optimization technique of 
Chan [IT] , to obtain an algorithm for the optimization problem, which runs in 0(n 3 log 5 n) ex- 
pected time. Our application of Chan's technique, described next, is somewhat non-standard, 
because each recursive step has also to handle global data, which it inherits from its ances- 
tors. 

Chan's technique, in its "purely recursive" form, takes an optimization problem that has 
to compute an optimum value w(P) on an input set P. The technique replaces P by several 
subsets Pi,...,P s , such that w(P) = mm{w(Pi), . . . , w(P s )}, and |Pj| < a\P\ for each % 
(here a < 1 and s are constants). It then processes the subproblems Pi in a random order, 
and computes minw(Pj) by comparing each w(Pi) to the minimum w collected so far, and 

by replacing w by w(Pi) if the latter is smaller jf| Comparisons are performed by the decision 
procedure, and updates of w are computed recursively. The crux of this technique is that 
the expected number of recursive calls (in a single recursive step) is only O(logs), and this 
(combined with some additional enhancements, which we omit here) suffices to make the 
expected cost of the whole procedure asymptotically the same as the cost of the decision 
procedure, for any values of s and a. Technically, if the cost D{n) of the decision procedure 

3 So the value of w keeps shrinking. 
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is J7(n 7 ), where 7 is some fixed positive constant, the expected running time is 0(D(n)) 
provided that 

(lns + l)a 7 <l. (1) 

However, even when (CQ) does not hold "as is" , Chan's technique enforces it by compressing 
I levels of the recursion into a single level, for I sufficiently large, so its expected cost is still 
0(D(n)). See [H] for details. 

To apply Chan's technique to our decision procedure, we pass to the dual space, where 
each point p G P is mapped to a plane p*, as done in the decision procedure. We obtain 
the set P* = {p* I p G P} of dual planes, and we consider its arrangement A = A(P*), 
where each cell r in A represents an equivalence class of planes in the original space, which 
separate P into the same two subsets of points P+, P~ . 

To decompose the optimization problem into subproblems, as required by Chan's tech- 
nique, we construct a (1/ ' g)- cutting of the dual space. We recall that, given a collection 
H of n hyperplanes in M. d and a parameter 1 < g < n, a (1 / g)- cutting of A(H) of size q 
is a partition of space into q (possibly unbounded) openly disjoint d- dimensional simplices 
Ai, . . . ,A q , such that the interior of each simplex A, is intersected by at most n/g of the 
hyperplanes of H. See [25] for more details. We use the following well known result [T3| [T^]: 

Lemma 3.2. Given a set H of n hyperplanes in M. d , a (1/ g)- cutting of A(H) of size 0(g d ) 
can be constructed in time 0{ng d ~ 1 ), for any g <n. 

Returning to our setup, we construct a (l/f?)-cutting for A(P*), for a specific constant 
value of g, that we will fix later, and obtain 0(g 3 ) simplices, such that the interior of each 
of them is intersected by at most n/g planes of P*. Each simplex Aj corresponds to one 
subproblem and contains some (possibly only portions of) cells 7~i, . . . ,Tfc of the arrangement 
A. We recall that each cell Tj represents an equivalence class of planes which separate P 
into two subsets of points P+ and P~. Hence, Aj represents a collection of such equivalence 
classes. All these subproblems have in common the sets (P*)^., (P*)^., consisting, respec- 
tively, of all the planes that pass fully above Aj and those that pass fully below Aj. (These 
sets are dual to respective subsets P£., of P, where Pj. is contained in all the sets P+, 
for the cells Tj, that meet Aj, and symmetrically for P^..) Note that most of the dual planes 
belong to (P*)J, U (P*)a • > the "undecided" planes are those that cross the interior of Aj, 
and their number is at most n/g. We denote the set of these planes as (P*)%. (and the set 
of their primal points as P£.)- 

To apply Chan's technique, we construct two segment trees on the arrangement of (P*)% , 
as described in Section I3TT1 Consider one of these segment trees, T + , that maintains the set 
of balls B + = {B r (p) I p G P/^}- Each cell Tj in Aj is represented by a leaf of T + . Each ball 
is represented as a collection of disjoint life-spans, with respect to a fixed tour of the cells of 
A((P*)° A ), which are stored as segments in T + , as described earlier. In addition, we compute 
the intersection of the balls centered at the points of P^., in 0(n log n) time, and store it at 
the root of T + . Note that, as we go down the recursion, we keep adding planes to (P*)^., 
that is, points to P^., and the actual set Pj. of points dual to the planes above the current 
Aj is the union of logarithmically many subsets, each obtained at one of the ancestor levels 
of the recursion, including the current step. However, we cannot inherit the precomputed 
intersections of the balls in these subsets of Pj. from the previous levels, since, as we go down 
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the recursion, Chan's technique keeps 'shrinking' the radius of the balls. Hence, each time 
we have to solve a decision subproblem, we compute the intersection of the balls centered 
at the points of Pj (collected over all the higher levels of the recursion) from scratch. (See 
below for details on the additional cost incurred by this step.) We build a second segment 
tree T~ that maintains the balls of B~ = {B r (p) \ p G P^}, in a fully analogous manner. 
The running time so far (of the decision procedure) is 0(nlogn + m 3 log 2 m), where m is 
the number of planes in (P*)% and n is the size of the initial input set P. 

To solve the decision procedure for a given subproblem associated with a simplex A,, we 
test, by going over all the root-to-leaf paths in T + and T~, whether there exists a cell r 
(overlapping Aj), for which the intersections of the spherical polytopes on the two respective 
paths in T + and T~ are nonempty (and, if nonempty, whether they both have nonempty 
interiors). The overall cost of this step, iterating over the 0(m 3 ) cells of A((P*)° A .) and 
applying the procedure from Section [37T1 for intersecting spherical polytopes, is 0(m 3 log 5 n). 

When the recursion bottoms out, we have two subsets P A ., and P A of 0(n) points, and a 
constant number of points in P A . Hence, we try the constant number of possible separations 
of P A . into an ordered pair of subsets Pi and P2, and, for each of these separations, we 
compute the two smallest enclosing balls of the sets Pj. U P\ and P A U P2 in linear time. 
If both P A . U Pi and P A U P2 can be covered by balls of radius r, for at least one of the 
possible separations of P£. into two subsets, then we have found a solution for the 2-center 
problem. (Discriminating between r* = r or r* < r is done as in Section I3~T1 ) 

We now apply Chan's technique to this decision procedure. Note that this application 
is not standard because the recursive subproblems are not "pure", as they also involve the 
"global" parameter n. We therefore need to exercise some care in the analysis of the expected 
performance of the technique. 

Specifically, denote by T(m, n) an upper bound on the expected running time of the 
algorithm, for preprocessing a recursive subproblem involving m points, where the initial 
input consists of n points. Then T(m, n) satisfies the following recurrence. 

~, v f \n(cg 3 )T(m/ g, n) + 0(m 3 log 5 n + nlogn), for m > g, , > 

^ \ 0(n), for m < g, ^ 

where c is an appropriate absolute constant (so that eg 3 bounds the number of cells of 
the cutting), and g is chosen to be a sufficiently large constant so that (JT|) holds (with 
s = eg 3 , a = 1/g, and 7 = 3). It is fairly routine (and we omit the details) to show that 
the recurrence ([2]) yields the overall bound 0(n 3 log 5 n) on the expected cost of the initial 
problem; i.e., T(n,n) = 0(n 3 log 5 n). We thus obtain the following intermediate result. 

Theorem 3.3. Let P be a set of n points in M 3 . A 2-center for P can be computed in 
0(n 3 log 5 n) randomized expected time. 

4 An Improved Algorithm 

4.1 An improved decision procedure T 

Consider the decision problem, where we are given a radius r and a parameter (3 > 0, and 
have to determine whether P can be covered by two balls of radius r, such that the distance 
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Figure 1: The points qi,q[,q2,q'2 prevent \c\C2\ from getting smaller. 



between their centers c\,C2 is at least j3r. (Details about supplying a good lower bound for 
/3 will be given in Section 14.21 ) By this we mean that there is no placement of two balls of 
radius r, which cover P, such that the distance between their centers is smaller than (3r; see 
Figure Q] 

This assumption is easily seen to imply the following property: Let C\ 2 denote the 
intersection circle of dB\ and dB 2 (assuming that B\ D B 2 7^ 0). Then any hemisphere 
v of <9£>i, such that (a) the plane n through c\ delimiting v is disjoint from Ci 2 , and (b) 
v and C\2 lie on different sides of n, must contain a point q of P, for otherwise we could 
have brought B\ and B 2 closer together by moving c\ in the normal direction of ir, into the 
halfspace containing c 2 (and Ci 2 ). See Figure [21 



Figure 2: The plane ir passes through c\ and is disjoint from C\ 2 - The hemisphere v delimited by tt, 
which lies on the side of ir not containing C± 2 , must contain a point q of P. 
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Figure 3: v\ is the leftmost point of the intersection circle C\ 2 . 



Guessing orientations and separating planes. We choose a set D of canonical orien- 
tations, so that the maximum angular deviation of any direction u from its closest direction 
in D is an appropriate multiple a of /3. The connection between a and (3 is given by the 
following reasoning. Fix a direction v G D so that the angle between the orientation of C\C2 
and v is at most a. Rotate the coordinate frame so that v becomes the x-axis. As above, 
let C\2 denote the intersection circle of dBi and dB 2 (assuming that the balls intersect). 
Let V\ be the leftmost point of C\ 2 (in the ^-direction); see Figure [3J If B\ and B 2 are 
disjoint (which only happens when \c\C2\ > 2r) we define V\ to be the leftmost point of B 2 . 
To determine the value of a, we note that (in complete analogy with Sharir's algorithm in 
the plane [31]) our procedure will try to find a yz-parallel plane, which separates c\ from 
v\. For this, we want to ensure that x[v\) — x(c\) > fir/A, say, to leave enough room for 
guessing such a separating plane. Let 9 denote the angle <V\C\C2 (see Figure HI). Using the 
triangle inequality on angles, the angle between c\v\ and the x-axis is at most 8 + a, so 
x{v\) — x(ci) > rcos(0 + a). Hence, to ensure the above separation, we need to choose a, 
such that cos(6 l + a) > (3/4. Since \cic 2 \ > fir, we have cos6 > (3/2. Hence, it suffices to 
choose a, such that 

a < cos cos — = sm sm — = B(p). 

4 2 2 4 y> ' 

With this constraint on a, the size of D is (1/a 2 ) = O (l/(3 2 ). 




x(ci) 



Figure 4: x(v±) — x{c\) > rcos(6 + a). 



We draw 0(1/(3) ?/2-parallel planes, with horizontal separation of /3r/4, starting at the 
leftmost point of P (with respect to the guessed orientation). One of these planes will 
separate v\ from c\. Thus, the total number of guesses that we make (an orientation in D 
and a separating plane) is 0(1/ (3 3 ). The following description pertains to a correct guess, in 
which the properties that we require are satisfied. (If all guesses fail, the decision procedure 
has a negative answer.) 

Reducing to a 2-dimensional search. By the property noted above, the left hemisphere 
v\ Q of dBi, delimited by the yz-parallel plane A through ci, must pass through at least one 
point q of P (see Figure [5]). 

Let Pl denote the subset of points of P lying to the left of A. Then P L must be fully 
contained in Bi and contain q. We compute the intersection K(P L ) = f]{B r (p) \ p G Pl} 
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Ac A 




Figure 5: The separating plane A and its parallel copy Ao through c±. The hemisphere v\ of dB\ to 
the left of Aq must contain a point q of P. 



in 0(n log n) time [9]. If K(Pl) is empty, then Pi cannot be covered by a ball of radius r 
and we determine that the currently assumed configuration does not yield a positive solution 
for the decision problem. Otherwise, since Pl Q #i, c\ must lie in K(Pl). Moreover, since 
q £ Pl lies on the left portion of dB\ , c\ must lie on the right portion of the boundary of 
K(Pl). Finally, since c\ lies to the left of A, only the portion ol of the right part of 8K(Pl) 
to the left of A has to be considered. If K(Pl) is disjoint from A then <jl is just the right 
portion of dK(P L ). Otherwise, o~l has a "hole", bounded by dK(P L ) R A, which is a convex 
piecewise-circular curve, being the boundary of the intersection of the disks B r (p) D A, for 
V e Pl- 

We partition o~l into quadratically many cells, such that if we place the center c\ of the left 
solution ball B>i in a cell r, then, no matter where we place it within r, B>i will cover the same 
subset of points from P. To construct this partition, we intersect, for each p £ P R = P\P L , 
the sphere dB r {p) with o L and obtain a curve 7 P on a L ; this curve bounds the portion of the 
unique face of 8K(Pl U {p}) within ctl- Hence, within K(Pl), it is a closed connected curve 
(it may be disconnected within o^, though). Let M denote the arrangement formed on o~l 
by the curves 7 P , for p £ Pr, and by the arcs of o~l- Apriori, M might have cubic complexity, 
if many of the 0(n 2 ) pairs of curves j a , 75, for a, b £ Pr, traverse a linear number of common 
faces of 0"^, and intersect each other on many of these faces, in an overall linear number of 
points. Equivalently, the "danger" is that the intersection circle C a b of a corresponding pair 
of spheres dB r (a), dB r (b), for a, b £ Pr, could intersect a linear number of faces of <jl (and 
each of these intersections is also an intersection point of 7 a and 7?,). See Figure El 

Complexity of M . Fortunately, in the assumed configuration, this cubic behavior is im- 
possible — C a b can meet only a constant number of faces of o^. Consequently, the overall 
complexity of M is only quadratic. This crucial claim follows from the observation that, for 
C a b to intersect many faces of <tl, it must have many short arcs, each delimited by two points 
on <jl and lying outside K(Pl). The main geometric insight, which rules out this possibility, 
and leads to our improved algorithm, is given in the following lemma. 

Lemma 4.1. Let X be a yz-parallel plane, which separates vi from c\. Let Pi C P be the 
subset of points of P to the left of X, and let Pr = P\Pl. Let C a b denote the intersection 
circle of dB r (a),dB r (b), for some pair of points a, b £ Pr, and let q £ Pl- If the arc 
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Figure 6: In a general setup (different than ours), an intersection circle of two balls (the dotted circle) 
may intersect a linear number of faces of 8K{Pl)- 



u = C a b \ B r (q) is smaller than a semicircle of C ao , then at least one of its endpoints must 
lie to the right of A. 

Proof. The situation and its analysis are depicted in Figure El To slightly simplify the 
analysis, and without loss of generality, assume that r = 1. Let h be the plane passing 
through a, b and q. Let c ao denote the midpoint of ab, and let w denote the center of the 
circumscribing circle Q of Aqab. Denote the distance \ab\ by 2x, and the radius of Q by y 
(so \wpi\ = \wp2\ = \wq\ = y). Note that c a b and w lie in h and that y > x. Observe that 
c a b is the center of the intersection circle C a b of dB r (a) and dB r {b). See Figure Eta). 

The intersection points z, z' of C a b arid dB r (q) are the intersection points of the three 
spheres dB r (a), dB r (b), and dB r (q). They lie on the line I passing through w and orthogonal 
to h, at equal distances a/1 — y 2 from w. See Figure Etb). (If y > 1 then z and z' do not 
exist, in which case C a b does not intersect dB r (q); in what follows we assume that y < 1.) 
Hence, within C a b, zz' is a chord of length 2a/1 — y 2 . In the assumed setup, z and z' delimit 
a short arc co of C a b, which lies outside B r (q), so points on the arc are (equally) closer to a 
and b than to q. 

Hence, the projection of the arc co onto h is a small interval ww', which lies on the bisector 
of ab in the direction that gets away from q; that is, it lies on the Voronoi edge of ab in the 
diagram Vor({a, b, q}) within h. See Figure Etc). Moreover, c a b also lies on the bisector, but 
it has to lie on the other side of w, or else the smaller arc co would have to lie inside B r (q). 
That is, c a b has to be closer to q than to a and b. Since A separates a and b from q, it also 
separates c a b from q. Moreover, the preceding arguments are easily seen to imply that wq 
crosses ab (as in Figure Eta)), which implies that A also separates q and w, so w has to lie 
to the right of A. Since z and z' lie on two sides of w on the line £, at least one of them has 
to lie on the same side of A as to (i.e., to the right of A). This completes the proof. □ 

Let a, b e Pr and consider those arcs of C a b which lie outside K(Pl) but their endpoints 
lie on <j£. Clearly, all these arcs are pairwise disjoint. At most one such arc can be larger 
than a semicircle. Let u be an arc of this kind which is smaller than a semicircle, and let 
q G Pl be such that one endpoint of u lies on dB r (q). Then lo' = C a b \ B r (q) is contained in 
co and therefore is also smaller than a semicircle. By Lemma 14.11 exactly one endpoint of uo' 
lies to the right of A (the other endpoint lies on cr^). Note that C a b cannot have more than 
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Figure 7: The setup in Lemma |4~T1 (a) the setup within the plane h; (b) the setup within C a b\ (c) 
ww' lies on the bisector of ab in the direction that gets away from q. 
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two such short arcs lying outside K(P L ), since, due to the convexity of C ao , only two arcs of 
C ao can have their two endpoints lying on opposite sides of A. Hence the number of arcs of 
C a b under consideration is at most 3, implying that 7 a and % intersect at most three times, 
and thus the complexity of M is 0(n 2 ), as asserted. 

Constructing and searching M. The next step of the algorithm is to compute M. We 
have already constructed 8K(Pl), in O(nlogn) time, and, in additional linear time, we can 
compute its portion ol to the left of A (we omit the straightforward details). We compute 
the intersection curve 7 P of B r (p) and <jl, for each p G Pr, in 0{n log n) time, by computing 
the intersection K(Pl U {p}), and obtaining the curve which bounds the portion of the 
unique face of dK(P L U {p}) within ol- If necessary, we also split 7 P into portions, such that 
each portion is contained in a different face of <jl- The total cost of computing all curves 
{7 P | p G Pr}, and spreading them along the faces of Ol, is 0(n 2 \ogn). Then, for each 
face / of a Li we consider the portions of all the arcs j p , for p e Pr, within /, and compute 
their arrangement (which is the portion of M which lies in /). To this end, we use standard 
line-sweeping [Tj5], to report all the intersections of n curves in the plane in 0((n + k) logn) 
time, where k = kf is the complexity of the resulting arrangement on /. Hence, the total 
cost of computing the portion of M on all the faces of ol is Sfeo- L ^(( n + kf)\ogn) = 
0(n 2 \ogn) + O(logn) • Ylfea L = 0{n 2 \ogn), since the complexity of M is 0{n 2 ). 

We next perform a traversal of the cells of M in a manner similar to the one used in 
Section [31 via a tour, which proceeds from each visited cell to an adjacent one. For each cell 
r that we visit, we place the center c\ of B\ in r, and maintain dynamically the subset of 
points of P not covered by B\. (Here, unlike the algorithm of Section [31 the complementary 
set P~ is automatically covered by B\ and there is no need to test it.) As before, when we 
move from one cell r to an adjacent cell t±, gains one point or loses one point. This 
implies that this tour generates only 0(n 2 ) connected life-spans of the points of P, where 
a life-span of a point p is a maximal connected interval of the tour, in which p belongs to 
P+. We can thus use a segment tree T M to store these life-spans, as before. Each leaf u of 
Tjvf represents a cell r of M, and the balls not containing r are those with life-spans that 
are stored at the nodes on the path from the root to u. Since M has a quadratic number of 
cells, T M has a total of 0{n 2 ) leaves. Arguing exactly as in Section [3TTI we can compute T M 
in overall 0(n 2 log 2 n) time, and the total storage used by Tm is 0(n 2 logn). 

As in Section [3TT] we next test, for each leaf u of Tm, whether the spherical polytopes 
along the path from the root to u have non-empty intersection. We do this using the 
parametric search technique described in Proposition 13. 1\ which takes 0(log 5 n) time for 
each path, for a total of 0(n 2 log 5 n). More precisely, as above, we also need to distinguish 
between r = r* and r > r*. We therefore stop only when both the intersection along the 
path and the cell of <jl corresponding to u are non-degenerate, and then report that r* < r. 
Otherwise, we continue running the above procedure over all paths of Tm, and repeat it for 
each of the 0(l//3 3 ) combinations of an orientation v and a separating plane A. If we find at 
least one (degenerate^) solution, we report that r* = r, and otherwise conclude that r* > r. 
Hence, the cost of handling Case 2, and thus also the overall cost of the decision procedure, 

4 Note that f]{B r (p) | p 6 P~} is non-degenerate if r is a 2-face or an edge. If r is a vertex we test for 
degeneracy as in the procedure in Section ETT1 Determining whether f\{B r {p) \ p 6 F+} is degenerate is also 
performed using that procedure. 
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is 0((l//3 3 )n 2 log 5 n). 



4.2 Solving the optimization problem 

We now combine the decision procedure T described in Section 14.11 with the randomized 
optimization technique of Chan [H] (as briefly described in Section 13721) . to obtain a solution 
for the optimization problem. 

The decision procedure T, on a specified radius r, relies on an apriori knowledge of a 
lower bound for the separation ratio \c\Ci\jr. To supply such a (3, let r denote the radius 
of the smallest enclosing ball of P, and observe that if there exist two balls Bi, B 2 of radius 
r covering P then the smallest ball B* enclosing B\ U B 2 must be at least as large as the 
smallest enclosing ball of P, so its radius must be at least r . Since this radius is (1 + /3/2)r 
(see Figured]), we have (1 + 0/2)r > r or > 2(r /r — 1). It follows that the running time 
of the decision procedure T is 




B* 



Figure 8: The smallest enclosing ball B* of Bi U In- 



dian's technique starts with a very big r (for all practical purposes we can start with 
r = r ) and shrinks it as it iterates over the subproblems. Therefore, running Chan's 
technique in a straightforward manner, starting with r = r , will make it potentially very 
inefficient, because the initial executions of T, when r is still close to r , may be too expensive 
due to the large constant of proportionality (not to mention the run at r itself, which the 
algorithm cannot handle at all). We need to fine-tune Chan's technique, to ensure that we do 
not consider values of r which are too close to r . To do so, we consider the interval (0,r ) 
which contains r*, and run an "exponential search" through it, calling T with the values 
Ti = ro (1 — 1/2 1 ), for % — 1, 2, . . ., in order, until the first time we reach a value r' = r, > r*. 
Note that 1 — r'/ro = 1/2 1 and 1/2* < 1 — r* /ro < 1/2*" 1 , so our lower bound estimates for 
the separation ratio /3 at r' and at r* differ by at most a factor of 2, so the cost of running T 
at r' is asymptotically the same as at r*. Moreover, since the (constants of proportionality in 
the) running time bounds on the executions of T at r\, . . . , form a geometric sequence, the 
overall cost of the exponential search is also asymptotically the same as the cost of running 
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r at r*. We then run Chan's technique, with r' as the initial minimum radius obtained so 
far. Hence, from now on, each call to T made by Chan's technique will cost asymptotically 
no more than the cost of calling T with r' (which is asymptotically the same as calling V 
with r*). 

Combining Chan's technique with the decision procedure V. To apply Chan's 
technique with our decision procedure, we use the same cutting-based decomposition as in 
Section [3T2l That is, we replace each point p G P by its dual plane p* G P*, and construct 
a (l/^)-cutting of A(P*), for some sufficiently large constant parameter g > 0. We then 
apply Chan's technique to the resulting subproblems (where each subproblem corresponds 
to a simplex Aj of the cutting), using the improved decision procedure T on each of them, 
and recursing into some of them, as required by the technique. As in Section [31 the recursion 
and the application of the decision procedure are not "pure" , because they need to consider 
also those planes that miss the current simplex. (Note that in the problem decomposition we 
use, for simplicity, the full 3-dimensional arrangement A(P*), of cubic size. This, however, 
does not affect the asymptotic running time, because we have only a constant number of 
subproblems, and Chan's technique recurses into only an expected logarithmic number of 
them.) Given a radius r, we compute the lower bound (3 = 2 — l) for the separation ratio 
J^i^i ; where c\, c 2 are the centers of the two covering balls, as above. Consider the application 
of r to a subproblem represented by a simplex Aj of the cutting. The presence of "global" 
points (those dual to planes passing above or below Aj) forces us, as in Section [372l to modify 
the "pure" version of Y described above. We use the same notations as in Section [3] 




A 7T 



Figure 9: h\ does not contain any point of P£.. 



We again rotate the coordinate axes, in 0(1/ P 2 ) ways (in the same manner as in the 
"pure" decision procedure), and draw 0(1/(3) yz-parallel planes, such that, at the correct 
orientation, one of these planes, A, separates C\ from v\ (if there is a solution for r). As in 
the pure case, we may assume that the x-span of P is at most 5r; a larger span is handled 
earlier. We assume, without loss of generality, that P^ C P>i, and that P£. C B 2 . Recall 
also that the points in the left halfspace h\ bounded by A are all contained in B\. Moreover, 
the plane tt containing the intersection circle Cyi is dual to a point it*, which has to separate 
(P*)X- from • Hence, all the points of P£. have to lie on the other side of tt, and in 

B> 2 , which is easily seen to imply that none of them can lie in h\. See Figure M We thus 
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verify that Pj. fl h\ = 0, aborting otherwise the guess of A. (Note that, in contrast, points 
of P^ can also lie to the right of A.) 

We now have a subset Pl C P^. of O(m) points to the left of A, which are assumed, 
together with the points of P^ , to be contained in B\. Note however that, for Lemma H~T1 to 
hold, we have to define 0£ only in terms of the points to the left of A. Therefore, we compute 
the surface a' L = dK(P L U (P^. D h\)) D h\ and search on it for a placement of the center 
Ci of B\. However, since the remaining points of P^ are also assumed to belong to Bi, we 
need to consider only the portion of a' L inside f){B r (p) \ p £ P^. \ h\}. Let o" L denote this 
portion. It is easy to compute a" L in 0(n\ogn) time. It is easily checked that c\ must lie on 
a'[ (if there is a solution for the current situation). So far, the cost of the decision procedure 
also depends (cheaply — see below) on the initial input size n, but the saving in this setup 
comes from the fact that it suffices to intersect the 0{m) spheres dB r (p), for p £ P^. \ h\, 
with a" L to obtain the map M, since only the points of P^. are "undecided". (The points of 
P£. are always placed in B 2 as already discussed.) 

Note that a'[ need not to be connected, so it may seem impossible to visit all the cells of 
M in a single connected tour. Nevertheless, we will be able to do it, in a manner detailed 
below. We thus build a segment tree T M to maintain the subset P'(ci) of points of P not 
covered by B\ . We build and query Tm as is done in Section 13.11 except for the following 
modifications. First, note that the points of P^ are assumed to be contained in B 2 . Thus, 
the points of P£ , that in the decision procedure were considered in building M, do not 
need to be considered as part of M now, rather it is enough to build the spherical polytope 
f]{B r (p) | p £ P^.} and place it at the root of T M . Second, we claim that M is of complexity 
0(mn). To see this, let C° denote the set of curves {dB r (p) flaj | p £ P^ }• Each pair 
of curves of C° can intersect each other in only a constant number of points, as proved in 
Section H~T1 Hence, the complexity of the arrangement of the 0(m) curves in C°, formed on 
a'[, is 0(m 2 ). However, a'[ itself is of complexity 0(n), and each edge of a'[ may intersect 
the curves of C° at 0(m) points. Hence, the complexity of the map M is 0(mn), but the 
number of its vertices that lie in the interior of the faces of M is only 0(m 2 ). 

To overcome the possible disconnectedness of a'[, we proceed as follows. We consider the 
(connected) network of the 0{n) edges of a' L , and intersect each of these edges with the m 
balls B r (p), for p £ P£.. We construct a tour of this network, which visits 0(mn) arcs along 
the edges of a' L , and append to this "master tour" separate tours of each face of a" L . We 
get in this way a single grand tour of the cells of M (which also traverses some superfluous 
arcs of a' L \ a'[), of length 0(mn), which has the incremental property that we need: Moving 
from any cell or arc of the tour to a neighbor cell or arc incurs an insertion or a deletion of 
a single point into/from P'(ci). 

Running time. For each cell of M we run the procedure described in Proposition 13.11 for 
determining whether the intersection of the corresponding spherical polytopes is nonempty 
(and whether it has nonempty interior). Therefore, solving each subproblem requires 0(mn log 5 n) 
time. The 0(mn log n) time required to build M, and the 0(n log ri) time required to con- 
struct the intersection of the balls in {B r (p) | p £ P£.}, are all subsumed in that cost. 
Repeating this for each of the 0(1/ f3 3 ) guesses of an orientation and a separating plane, 
results in O ((l//3 3 )mn log 5 nj rnning time. When the recursion bottoms out, we handle it 
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the same way as in Section 13721 

Arguing similarly to the less efficient solution, we obtain the following recurrence for 
the maximum expected cost T(m, n) of solving a recursive subproblem involving m "local" 
points, where n is the number of initial input points in P. 

Tim ) < I m ( c £ 3 ) r ( m M n ) + ((1 /(3 3 )mn log 5 n) , for m > q, , , 

^ \ 0(n), for m < g, ^ ' 

where c is an appropriate absolute constant (as in Section [372]), q is the parameter of the 
cutting, chosen to be a sufficiently large constant (to satisfy (JT]), as above, with 7 = 2), 
and P = 2 (r /r' — 1), where r' is the value of r at which the initial exponential search is 
terminated. 

It can be shown rather easily (and we omit the details, as we did in the preceding section), 
that the recurrence ([3]) yields the overall bound O ((l//3 3 )n 2 log 5 n) on the expected cost of 
the initial problem; i.e., 

T(n,n) = O ((l//3 3 )n 2 log 5 n) . 
We thus finally obtain our main result: 

Theorem 4.2. Let P be a set of n points in M 3 . A 2-center for P can be computed in 
0((n 2 log 5 n) / (1 — r*/r ) 3 ) randomized expected time, where r* is the radius of the balls of 
the 2-center for P and r is the radius of the smallest enclosing ball of P. 



5 Efficient Emptiness Detection of Intersection of Spher- 
ical Polytopes 

In this section we describe an efficient procedure for testing emptiness (and non-degeneracy) 
of the intersection of spherical polytopes, as prescribed in Proposition 13.11 Let S be a 
collection of spherical polytopes, each defined as the intersection of at most n balls of a fixed 
radius r. Fix a spherical polytope S G S. To simplify the forthcoming analysis, we assume 
that the centers of the balls involved in the polytopes of S are in general position, meaning 
that no five of them are co-spherical, and that there exists at most one quadruple of centers 
lying on a common sphere of radius r. As is well known, each ball b participating in the 
intersection S contributes at most one (connected) face to dS (see [32]). The vertices and 
edges of S are the intersections of two or three bounding spheres, respectively (at most one 
vertex might be incident to four spheres). Hence dS is a planar (or, rather, spherical) map 
with at most \S\ faces, which implies that the complexity of dS is 0(|5|). 

We preprocess S into a point-location structure. We first partition dS into its upper 
portion dS + and lower portion dS~. We project vertically each of dS + and dS~ onto the 
xy-plane and obtain two respective planar maps M + and M~ (see Figure [TU]) . For each face 
C of each map we store the ball b that created it; that is, £ is the projection of the (unique) 
face of dS that lies on db. The rry-projection S* of S is equal to both projections of dS + , 
dS~ , and is bounded by a convex curve E* that is the concatenation of the x|/-projections 
of certain edges of S and of portions of horizontal equators of some of its balls. 

We apply the standard point-location algorithm of Sarnak and Tarjan [33] to each of 
the maps M + ,M~. That is, we divide each planar map into slabs by parallel lines (to the 
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Figure 10: Projecting dSi vertically onto the xy-plane (left), and the point location structure for the 
resulting map Mf (right). 



y-axis) through each of the endpoints (and locally x-extremal points) of the arcs obtained 
by projecting the edges of OS, including the new equatorial arcs. Using the persistent search 
structure of [33], the total storage is linear in \S\ and the preprocessing cost is 0(\S\ log \ S\), 
where \S\ is the number of balls forming S. To locate a point go i n M + (or in M~), we first 
find the slab in the x-structure that contains go, and then find the two curves between which 
go lies in the y-structure|^ 

To determine whether g G S*, we locate the face ( + (resp., £~) of the map M + (resp., 
M~) that contains g, as just described. Each of these faces can be a 2-face, an edge or a 
vertex. We therefore retrieve a set B + (resp., B~) of the one, two, or three or four balls 
associated (respectively) with the 2-face, edge or vertex containing q. (We omit here the 
easy construction of witness balls when the faces ( + and (~ are not associated with any ball, 
that is, q £ S*.) 

Let B denote the set B + U B~ . We observe that q G S* if and only if the z- vertical 
line X q through q intersects S. Moreover, we have, by construction, X q fl S = X q fl (f) B). 
Hence q G S* if and only if s := \ q fl (P| B) ^ 0. Clearly, if we put N = J2ses 1^1' then the 
preprocessing stage takes a total of 0(iV log n) time and requires O(N) storage. 

Next, let Si, . . . , S t be t < logn spherical polytopes of S, for which we want to determine 
whether K = Hi=i &i * s nonempty (and, if so, whether it has nonempty interior). We solve 
this problem by employing a technique similar to the multi-dimensional parametric searching 
technique of Matousek (26] (see also [H [101 [T71 Ej). We solve in succession the following 
three subproblems, n (g), where g is a point in the xy-p\ane, ITi (Z) , where I is a y-parallel 
line in the xy-plane, and n 2 , over the entire xy-plane. In the latter problem we wish to 
to determine whether the xy-projection K* of K is nonempty. During the execution of 
the algorithm for solving Il 2 , we call recursively the algorithm for solving IIi(Z), for certain 
^/-parallel lines I C 1R 2 , and we wish to determine whether K* meets I. If so, then II2 is 
solved directly (with a positive answer). Otherwise, we wish to determine which side of I, 
within R 2 , can meet K* (since K* is convex, there can exist at most one such side). The 
recursion bottoms out at certain points q G /, on which we run n (g) to determine whether 
K* contains g. If so, then H\(l) is solved directly (with a positive answer). Otherwise, we 
determine which side of g, within I, can meet K*, and continue the search accordingly. 

Our solutions to the subproblems Hfc, < k < 2, are based on generic simulations of the 

5 All these standard details are presented to make more precise the infrastructure used by the 
higher-dimensional routines III and II2. 
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standard point-location machinery of Sarnak and Tarjan [33] mentioned above. In each of 
the subproblems, if we find a point in / D K*, for the respective point, line, or the entire 
xy-plane /, we know that K 7^ and stop right away. If / D K* = 0, we want to "prove" 
it, by returning a small set of witness balls b\, . . . ,b y , where, for each j, bj is one of the 
balls that participates in some spherical polytope Si (so bj D Si), so that their intersection 
K = f)? =1 bj satisfies / fl Kq = (where, as above, Kq is the ^-projection of Kq). If 
Kq = then K = too and we stop. Otherwise (when / is a line or a point), Kq determines 
the side of / (within M 2 if / is a line, or within the containing line I if / is a point) that 
might meet K*\ the opposite side is asserted at this point to be disjoint from K*. We use 
this information to perform binary search (or, more precisely, parametric search) to locate 
K* within the flat, from which we have recursed into /. The execution of the algorithm 
for solving II2 will therefore either find a point in K or determine that K — 0, because 
it has collected a small (as we will show, polylogarithmic) number of witness balls, whose 
intersection, which has to contain K, is found to be empty. 

Solving Ho(q) for a point q. Here we have a point g e R 2 and we wish to determine 
whether q e K*. To do so, we locate q in each of the maps Mf (the xy- projection of dS^~) 
and M[~ (the ^-projection of dS~), for each % — 1, . . . , t. If q lies outside the projection of 
at least one polytope Si then q ^ K*, and we return the witness balls that prove that q ^ S*. 
Otherwise, as explained above, each point location returns a set £>j of 0(1) witness balls for 
S^ We compute the t line segments Sj = X q fl (f]Bi), for each i — 1, . . . ,t, where X q is, as 
above, the ^-vertical line through q. We then have Kq := X q fl K = Hi=i s ^ so ^ suffices to 
compute this intersection (in 0(t) time) and test whether it is nonempty. If Kq is nonempty, 
then we have found a point q' in K. Otherwise, we return the set Bq = [J{Bi | 1 < % < t} of 
up to 51ogn balls as witness balls for the higher-dimensional step (involving the y-parallel 
line containing q). 

The time complexity for solving Ilo(g) is 0(log n), since it takes O(logn) time to com- 
pute, for each of the O(logn) spherical polytopes Si, the intersection \ q fl Si. 

Solving IIi(Z) for a line I. Here we have a y-parallel line I C M 2 and we wish to determine 
whether K* meets I. We first locate I in each of the planar maps Mf and M~ of each Si, 
and find the slabs ^ and ip~ , which contain / (in some cases I is the common bounding 
line of two adjacent slabs ip[ and ip" of Mf or of , so we retrieve both slabs). We then 
run a binary search through the y-structure of each of the obtained slabs to find a point in 
K* n I, if one exist. In each step of the search, within some fixed slab ipo, we consider an 
arc 7 of the ^/-structure, and determine whether K* meets / above or below 7 (within IR 2 ), 
assuming K* fl I ^ 0. To this end, we find the intersection point go = I H 7, and run the 
algorithm for solving n (g ) (see Figure [TT]). If q e K* , then we have found a point q' in K, 
and we immediately stop. Otherwise, we have a set Bq of up to 5 logn balls returned by the 
algorithm for solving n (go)- We test whether the xy-projection Kq of H^o intersects I. If 
Kq fl I = 0, then (due to the convexity of K) we know which side of I (within IR 2 ) meets K* , 
and we return Bq as a set of witness balls for the higher- dimensional (planar) step. Otherwise 
(again due to the convexity of K), we know which side of 7, within I, meets K*, and we 
continue the search through the ^/-structure of ipo on this side. We continue the search in 
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this manner, until, for each Si, we obtain an interval £j of Z between two consecutive arcs of 
the ^-structure of ipo, which meets K* (assuming K* nZ 7^ 0). Let S denote the collection of 
all these intervals. Clearly, K* D Z C P)H. We find the lowest endpoint among the top 
endpoints of the intervals in S and the highest endpoint E + among the bottom endpoints 
of the intervals in 5, and test whether E~ is above E + . If so, we consider the set B\ of up 
to 10 log n witness balls returned by the algorithms for solving U (E~) and Uq(E + ). If the 
xy-projection K* of f] B\ intersects Z, then K* meets Z and we stop immediately, for we have 
found that K is nonempty. Otherwise, we know which side of Z (within IR 2 ) can meet K*, 
and we return B\ as a set of witness balls for the higher (planar) recursive level. If E~ is not 
above E + , then K* fl Z = and we return B\ as a set of witness balls for the higher (planar) 
recursive level as welljf] 



to 
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Figure 11: The line I on which we run Ili(Z). The point qo on which we run Ilo(go) is the intersection 
point of I with some arc 7. 

A naive implementation of the above procedure takes 0(log 4 n) time, since for each of 
the O(logn) spherical polytopes Si we run a binary search through the ^-structure of at 
most two slabs of each of the maps Mf and M~ , and in each of the binary search steps, we 
run the algorithm for solving Ho(go) for some point go- The other substeps take less time. 
However, we can improve the running time by implementing it in a parallel manner and 
simulating the parallel version sequentially with a smaller number of calls to n . 

We only parallelize the binary searches through the y-structure of each Mf and M[~ , 
since the other substeps take less time. To this end, we use O(logn) processors, one for 
each of the planar maps Mf and M, - , and we run in parallel the binary search through the 
^/-structure of each planar map using O(logn) parallel steps. In each parallel step we need 
to "compare" O(logn) arcs with K* (one arc for each of the planar maps M± , Mf). We 
therefore intersect each such arc with Z and obtain a set Q of O(logn) intersection points. We 
then run a binary search through the points of Q (to locate K*) using Ho- This determines 
the outcome of the comparisons of each of the arcs with K*, and the parallel execution 
can proceed to the next step. Applying this approach to each of the O(logn) parallel steps 
results in an 0(log 3 n log log n)-time algorithm for solving H\(l). However, we can slightly 
improve this bound further using a simple variant of Cole's technique [18]. More precisely, 

With some care, the number of witness balls can be significantly reduced. We do not go into this 
improvement, because handling the witness balls is an inexpensive step, whose cost is subsumed by the cost 
of the other steps of the algorithm. 
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in each parallel step we have a collection Q of O(logn) weighted points, one for each map, 
which we need to compare with K*. We select the (weighted) median point go of Q and run 
Ho (go). This determines the outcomes of the comparisons between K* and each of the points 
in Q which lie to the opposite side of go to the side containing K* . Points in Q which lie in 
the same side of go as K*, in level j of the parallel implementation, are given weight 1/4 J_1 
and we try to resolve their comparison to K* in the next step. An easy calculation (simpler 
than the one used by Cole) shows that this method adds only O(logn) steps to the O(logn) 
parallel steps of the searches, and now in each parallel step we perform only one call to n 
(see [18] for more details). Therefore, the total running time of IIi(Z) is 0(log 3 n). 

Solving Il 2 . We next consider the main problem Il 2 , where we want to determine whether 
K* 7^ (i.e., whether K ^ 0). We use parametric searching, in which we run the point 
location algorithm that we used for solving IIo, in the following generic manner. 

In the first stage of the generic point location, we run a binary search through the slabs 
of each of the planar maps Mf and M~ , for i — 1, . . . , t. In each step of the search through 
any of the maps, we take a line Iq delimiting two consecutive slabs of the map, and run the 
algorithm for solving IIi(Zo), thereby deciding on which side of Iq to continue the search. At 
the end of this stage, unless we have already found a point in K or determined that K is 
empty, we obtain a single slab in each map that contains K*. Let ip denote the intersection 
of these slabs, which must therefore contain K* (unless K is empty). The cost of this part 
of the procedure is 0(log 5 n). 

In the next stage of the generic point location, we consider each map Mf or M~ (for 
simplicity we refer to it just as Mi) separately, and run a binary search through the y- 
structure of its slab ipi that contains ip. In each step of the search we consider an arc 7 of 
the ^-structure, and determine which side of 7 (within the slab ift), can meet K*, assuming 
that ip fl K* 7^ 0; if 7 fl K* 7^ we will detect it and stop right away. Before describing in 
detail how to resolve each comparison with an arc 7, we note that this results in O(logn) 
comparisons of arcs 7 to K* for each of the O(logn) planar maps Mf and M~ . However, 
we can reduce the number of comparisons to O(logn) in total, by simulating (sequentially) 
a parallel implementation of this step, as follows. There are O(logn) parallel steps, and in 
each step we execute a single step of the binary search in each of the maps Mf , M~ . In each 
parallel step we need to compare K* to a set G of O(logn) arcs, one of each of the planar 
maps Mj + , M~ . Consider the portion A'(G) of the arrangement A{G) of the arcs in G which 
lies in t/j. Let L(G) denote the set of 0(log 2 n) y-parallel lines which pass through the vertices 
of A'(G). We run a binary search through the lines of L(G), using calls to the algorithm 
for rix to guide the search, to locate K* amid these lines, in a total of O (log 3 n log log n) 
running time. This step (if it did not find a line crossing K*) may trim ip to a narrower 
slab ip' in which K* must lie if K* ^ 0. Put G' = {7 fl tp' \ 7 G G}, and observe that 
the arcs of G' are pairwise disjoint and form a sorted sequence in the y-direction. We then 
perform a binary search through the arcs in G', using O(loglogn) comparisons to K* . Each 
comparison is carried out in 0(log 4 n) time, in a manner detailed below. Once the binary 
search is terminated, we can determine the outcomes of the comparisons between K* and 
each of the arcs in G' and proceed to the next parallel step. Applying this approach to 
each of the O(logn) parallel steps results in an 0(log 5 n log log n)-algorithm for solving II2. 
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We again use an appropriate variant of Cole's technique to improve the running time by a 
log log n factor, in a manner similar to the one described in the solution of III. 

To carry out a comparison between an arc 7 G G' and K*, we act under the assumption 
that 7 fl K* 7^ 0, and try to locate a point of 7 n K* in each of the other maps. Suppose, to 
simplify the description, that we managed to locate the entire 7 in a single face of each of 
the other maps Mj~, M~ . This yields a set B of 0(t) balls, so that a point v 6 7 lies in K* 
if and only if it lies in the xy-projection Kq of f]B. We then test whether 7 intersects Kq. 
If so, we have found a point in K and stop right away. Suppose then that Kq fl 7 = 0. If 
Kq dip' = then K must be empty, because we already know that K* C ip'. If Kq dip' 7^ 0, 
then we know on which side of 7 to continue the binary search in (the portion within ip' of) 




k 

Figure 12: Comparing 7 n K* with 5, The outcome of IIi(Zo) determines (a) the side of 5 in which 
the search in ipj should continue, and (b) the portion of 7 which can still meet K*. The subslab 1// is 
drawn shaded. 



In general, though, 7 might split between several cells of a map Mj, where Mj denotes, as 
above, one of the maps Mj~ or M~ . This forces us to narrow the search to a subarc of 7, in 
the following manner. We run a binary search through the ^/-structure of the corresponding 
slab ipj of Mj, which contains ip', and repeat it for each of the maps Mj. In each step of the 
search, we need to compare 7 (or, more precisely, some point in 7 fl K*) with some arc 5 
of ijjj, which we do as follows. If 7 lies, within ip', completely on one side of 5, we continue 
the binary search in ifjj on that side of 5. If 7 intersects 5, we pick an intersection point v 
of 7 and 5, pass a y-parallel line Iq C 1R 2 through v, and run the non-generic version of the 
algorithm to solve IIi(Zo)- (See Figure fT2l) As before, if Iq fl K* 7^ we detect this and stop. 
Otherwise, we know which of the two portions of 7, delimited by v, can intersect K*. We 
repeat this step for each of the at most four intersection points of 7 and 5 (observing that 
these are elliptic arcs), and obtain a connected portion 7' of 7, delimited by two consecutive 
intersection points, whose relative interior lies completely above or below 5, so that 7 fl K*, 
if nonempty, lies in 7'. This allows us to resolve the generic comparison with 5, and continue 
the binary search through ipj. (On the fly, each comparison with a line Iq narrows ip' still 
further.) 

To make this procedure more efficient, we perform the binary searches through the slabs 
ipj in parallel, as follows. As before, we run in parallel the binary searches through each of 
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the slabs ipj using O(logn) parallel steps. In each parallel step we need to compare a set 
D of O(logn) arcs to 7, one arc 5 from each planar map M y We intersect each of the arcs 
in D with 7 and obtain a set Z of O(logn) intersection points. Let Lz denote the set of 
the O(logn) y-parallel lines which pass through the points of Z. We run a binary search 
through the lines of Lz-, using calls to the algorithm for LTi to guide the search, in a total 
of O (log 3 n log log n) running time. We obtain a connected portion 7' of 7, delimited by two 
consecutive intersection points of Z , whose relative interior lies completely above or below 
each 5 G D, so that / -f(lK*, if nonempty, lies in 7'. This allows us to resolve each comparison 
between K* and an arc 5 G D, assuming that ^(IK* 7^ 0, and we continue the binary search 
through each Mj in the same manner. 

We again use a variant of Cole's technique [18] to slightly improve this bound further. 
In each parallel step we have a collection Z of 0(logn) weighted points, each of which is an 
intersection point of 7 with some arc S from one of the planar maps Mj, and we need to 
compare each of the points of Z with K* . Let D denote the set of these active arcs. 

Note that each arc 5 participating in this step contributes (at most) four points to Z, 
for a total of at most 4|D| points. We perform three steps of a (weighted) binary search on 
the points of Z, where each step takes the weighted median z of an appropriate portion 
of Z, and calls n 1 (/ ), where Iq is the vertical line through zq. These IL-steps resolve the 
comparisons with K* of all but 1/8 of the points of Z, that is, at most (1/8) • 4|D| = \D\/2 
points of Z are still unresolved. 

In other words, after the three calls to the algorithm for solving LTi (in the first parallel 
step of the execution), we can determine the outcomes of the comparisons of at least half of 
the arcs in D with K*. We can then proceed in this manner and apply Cole's technique (as 
before), by using only a constant number of calls to II 1 in each of the O(logn) parallel steps 
of searching in all the maps. This reduces a log log n factor from the bound of the running 
time, so it is only 0(log 5 ri) time. 

When these searches terminate, we end up with a 2-face in each Mj, in which 7 fl K* lies 
(if nonempty), and we reach the scenario described in a preceding paragraph. As explained 
there, we can now either determine that K ^ 0, or that K = 0, or else we know which side 
of 7, within ipi (or, rather, within ip') can contain K*, and we continue the binary search 
through ipi on that side. 

When the binary search through ipi terminates, we have a 2-face Q of Mj, where K* 
must lie, and we retrieve the ball bi corresponding to Q. We repeat this step to each of 
the maps Mf and M~ of each of the t spherical polytopes Si, and obtain a set B\ of 2t 
balls. In addition, the searches through the maps Mf and M~ may have trimmed ip' to a 
narrower strip ip", and have produced a set B2 of witness balls, so that the xy-projection of 
their intersection lies inside ip". B2 may consist of a total of 0(t 3 log 2 n) witness balls, as 
is easy to verify. In addition, the second-level searches produce an additional collection B' 2 , 
consisting of balls corresponding to faces of the maps Mj~ and M~ , in which the second-level 
searches have ended; their overall number is 0(t 2 logn). Put K2 = f](B\ U B2 U B' 2 ). Hence 
K ^ if and only if K 2 ± 0. 

As already noted, the overall running time of the emptiness detection is 0(log 5 n). 

So far, we have only determined whether K is empty or not. However, to enable the 
decision procedure to discriminate between the cases r* = r and r* < r we need to refine 
the algorithm, so that it can also determine whether K has nonempty interior (we refer to 
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an intersection K with this property as non- degenerate). To do so, we make the following 
modifications to the algorithm described above. Each step in the emptiness testing procedure 
which detects that K ^ obtains a specific point w that belongs to K. Moreover, w belongs 
to the intersection K\ of polylogarithmically many witness balls, and does not lie on the 
boundary of any other ball. This is because each of the procedures n ,rii, or n 2 locates 
the :q/-projection w* of w (which, for Hi and n 2 is a generic, unknown point in K) in each 
of the maps M^~, M~ , i = 1, . . . , t, and the collection of the witness balls gathered during 
the various steps of the searches contains all the balls that participate in the corresponding 
spherical polytopes Si on whose boundary w can lie. Thus, when we terminate with a point 
w G K, we find, among the polylogarithmically many witness balls, the at most four balls 
whose boundaries contain w (recall our general position assumption), and test whether their 
intersection is the singleton {w}. It is easily checked that this is equivalent to the condition 
that K is degenerate. 

This completes the description of the algorithm, and concludes the proof of Proposi- 
tion O 

6 Discussion and Open Problems. 

In this paper we presented two algorithms for computing the 2-center of a set of points in 
R 3 . The first algorithm takes near-cubic time, and the second one takes near-quadratic time 
provided that the two centers are not too close to each other. Note that our second algorithm 
may be slightly revised, so that it receives, in addition to P, a parameter e > as input, 
and returns a solution for the 2-center problem for P, if e < 1 — To this end, we run 
the exponential search until we reach a value of r with 1 — ^ < e. If along the search we 
have found a value of r such that r > r*, we stop the search and run Chan's technique with 
the constraint that r* < r, as above. Otherwise, we have r* > ro(l — e) and we may return 
the smallest enclosing ball of P as an e-approximate solution for the 2-center problem. This 
way, we ensure that the running time of our algorithm is 0(e~ 3 n 2 log 5 n). 

An obvious open problem is to design an algorithm for the 2-center problem that runs 
in near-quadratic time on all point sets in M 3 . Another interesting question is whether 
the 2-center problem in M 3 is 3SUM -hard (see [22] for details), which would suggest that a 
near-quadratic algorithm is (almost) the best possible for this problem. 
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